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Heating of figure-8 ions by odd-parity rotating mag- 
netic fields [RMFo] applied to an elongated field- 
reversed configuration (FRC) is investigated. The largest 
energy gain occurs at resonances (s = lor/lo) of the 
RMFo frequency, uir, with the figure-8 orbital frequency, 
w, and is proportional to for s — even resonances and 
to s for s — odd resonances. The threshold for the tran- 
sition from regular to stochastic orbits explains both the 
onset and saturation of heating. The FRC magnetic ge- 
ometry lowers the threshold for heating below that in the 
tokamak by an order of magnitude. 

Heating, i.e., stochastic energy gain, of charged parti- 
cles by time-varying fields is a complex and fundamen- 
tal phenomenon critically important to as diverse areas 
of plasma physics as fusion research [T] and plasma pro- 
cessing. Well known are the effects of simple resonances 
and particle collisionality on the heating of magnetized 
plasmas. Far less well explored is the role of an inho- 
mogeneous static magnetic-field geometry. Because of 
its relevance to space plasmas [3], plasma processing ^4j, 
and magnetic-confinement controlled- fusion researchjS], 
the field-reversed configuration (FRC, see Fig. 1) - with 
its poloidal field nulls, lack of toroidal field, and strong 
field gradients - is an important system in which to ex- 
plore the effects of magnetic field geometry on particle 
dynamics under the influence of time- varying fields. 

Even with axial symmetry, a static FRC allows 
charged-particle orbits that are regular or ergodic'llj. 
First studies of single-particle orbits in FRCs assumed 
time invariance and spatial symmetries that reduced the 
problem to one or two dimensions, allowing Kolmogorov- 
Arnold-Mosher (KAM) surfaces to exist [2 and limiting 
excursions in phase space. The addition of a rotating 
magnetic field (RMF)[TT] breaks the angular invariance 
of the FRC, creating a three-dimensional system without 
bounding KAM surfaces and opening the possibility for 
large excursions in phase space and energy. These excur- 
sions can have beneficial results, such as ion heatingPT], 
or detrimental ones, such as loss of confinement. In this 
paper we present studies of ion orbits in FRCs with RMF 
applied: the goal is to understand the threshold for chaos 
and the role of resonances in the non-linear growth and 
subsequent saturation of ion energy. We restrict atten- 



tion to the novel odd-parity RMFs [RMFo] because of 
field closure and encouraging recent experimental results. 
We show that the same mechanism is responsible for the 
initial ion heating and its ultimate saturation. 

Studies of stochastic ion heating by perpendicularly 
propagating electrostatic waves in tokamaks were per- 
formed with similar Hamiltonian techniques and research 
goals. The results we report are markedly different be- 
cause of fundamental differences in the magnetic field 
geometry of the two devices. 

Earlier papers U , which used the RMF numerical code 
to investigate RMFqS applied to FRCs, showed that 
the relevant frequency range for ion heating was broad, 
~ 0.2 — 2, where = w/j/Wci, is the RMFo fre- 
quency, bJci = qBa/mc is the ion-cyclotron frequency in 
the axial field at the FRCs center, Ba, m is the ion mass, 
and q is the ion charge. These papers reported signifi- 
cant ion heating even for low relative RMF amplitude, 
Br'. Bfi/Ba ~ 5 X 10"'*. Phase de-coherence of ion or- 
bits, with respect to the periodic electric fields created 
by the RMFo, is a necessary condition for ion heating. 
Strong gradients and regions of field reversal in the FRC 
provide locations for possible phase de-coherence. For a 
10-cm FRC having an ion density of 10*^ cm"'^ and an 
ion energy of 100 eV, Coulomb collisions will be 10 x less 
frequent than the stochastic effects described herein [7]. 

The question arose whether, in spite of the existence 
of strong field gradients, ion-cyclotron resonances (ICRs) 
were important to ion heating. We show that ICRs are 
important, but with significant differences from the stan- 
dard ICR picture. More rapid heating occurs at low 
Br/B^ for figure-8 orbits (see Fig. lb)) than for cy- 
clotron orbits, though the latter have a more clearly 
resonant interaction with RMFq. Figure-8 orbits cross 
the field-reversal and strong-gradient regions (twice) ev- 
ery orbit cycle, possibly losing phase coherence at each 
traversal. In contrast, cyclotron orbits may only incur 
phase de-coherence at the less frequent excursions to the 
axial extremes of their orbits. Betatron orbits have a less 
non-linear nature and hence are also less well heated than 
figure-8 orbits. Because figure-8 orbits are representative 
of a large fraction of ions in hot fusion FRC plasmas and 
because they represent the physically interesting situa- 
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tion of motion in a double potential well[Tn], we focus 
on them. The studies presented herein also clarify why 
high-energy orbits tend to interact regularly with RMFo, 
leading, importantly, to a saturation of ion heating by 
RMFo and a method for tuning ion energy. 

We follow Ref. [7] by using the same equations for the 
RMFo and a Solov'ev equilibrium for the FRC, with the 
notation, see Fig. la): R = FRC separatrix radius in 
midplane; Z = FRC axial half length; k = Z/R, FRC 
elongation; r = radial coordinate; z = axial coordinate; 
(j) = azimuthal coordinate; pi = canonical momenta; P = 
2p^/qBaR^ , normalized p^; A = vector potential of RMF 
and FRC; ip = (p — uJut; k = Itt/kR = axial wave number 
of the RMFo] I = axial mode number; and time, r, in 
units of 27r/wci- 

The shape of the effective potential-energy surface on 
which an ion moves depends on P and z[? ? ]: figure-8 
orbits may be confined to the z = subspace, or to a 
potential- well minimum above or below 2; = 0, or may 
oscillate across z = 0. Orbits confined to the z — 
subspace are amenable to an analytic analysis and are 
the appropriate choice to analyze because of the RMFo^s 
electric fields, £r and 5^, there. Each cross-section in z 
is either a double potential well, allowing both cyclotron 
and figure-8 orbits, or a raised potential, corresponding 
to betatron orbits. Cyclotron orbits feel a force towards 
larger \z\, thus eventually enter a region where the barrier 
between the double wells is low enough for them to tra- 
verse, thereby becoming figure-8 orbits. Since cyclotron 
orbits interact regularly with RMF, except at these ax- 
ial extremes, their random fluctuations in energy appear 
less frequently than for orbits which are always figure-8. 
Moreover, figure-8 orbits have greater radial excursions, 
hence gain more energy from the radial electric field of 
the RMFo- It follows that the heating of figure-8 orbits 
is an upper limit for the heating of all ions in the FRC 
and that the threshold for heating is highest in the z = 
subspace. Extensive numerical simulations, performed 
with the RMF code, confirmed this. 

We first examined whether the broad Q range for heat- 
ing is due to resonances at the fundamental ICR fre- 
quency. As described below, the answer is no. Instead 
high-harmonic resonances occur because the frequency 
of the figure-8 orbit is highly nonlinear. (High-harmonic 
resonances have recently been observed in an RMF ex- 
periment.) As the energy of a figure-8 orbit decreases, the 
ratio s = uir/cj increases because the ion's frequency, w, 
slows down as it gets closer to the phase-space separatrix 
created by the hump in the double potential well. A set 
of resonances with the RMFo occurs at integral values 
of s. 

Figure 2a) shows the i?MF-code-calculated time de- 
pendence of ion energy for two values of Bb, for a 1-keV 
ion initiated in a figure-8 orbit with zero axial velocity 
in the z = subspace of an FRC having R = 10 cm, 
K = 5, il = 0.9, and Ba = 20 kC Regular motion, with a 



clear s = 5 component, is seen for Bn = 2 G. The energy 
fluctuations are small (~ 15%) for Bn = 2 G and large, 
> 100%, for Br = 20 G. For Br = 2G, the fast Fourier 
transform (FFT) of ion energy. Figure 2 b), shows sharp 
peaks in frequency space, indicative of regular motion. 
The separation between peaks, A/, is 0.207 ± 0.002, in 
units of ujci- For Br — 20 G, the FFT shows broadband 
noise at a 30 x higher absolute level. Under these con- 
ditions, betatron (also shown in Fig. 2b)) and cyclotron 
orbits (not shown) display regular motion — sharp peaks 
in their FFTs — with energy fluctuations less than 7%, 
even for Br = 20 G. Figure 2c) shows A/ vs E, initial 
energy normalized to En = q^B'^R^ /2m, for three values 
of initial P and low Br. Below P = 0.25 orbits may be 
cyclotron or figure-8; above P — 0.25 orbits are betatron. 
A logarithmic drop in A/ is seen for P = 0.1 and 0.2 at 
an energy corresponding to the phase-space separatrix 
energy, at the transition of cyclotron into figure-8 orbits. 
Little change in A/ occurs for betatron orbits. 

In heating, the variance of energy, and therefore the 
maximum energy, E^ax, will increase with time. Figure 
3a) shows Emax attained by an initially figure-8 orbit 
as a function of time for four Br values and the same 
FRC parameters as in Fig. 2. E^ax displays saturation 
behavior quickly, implying phase coherence growing with 
increasing energy. The threshold for heating is at Br ^ 3 
G, above which E^nax grows a B]^^. Figure 3b) compares 
E.nax VS Br for figure-8 (P = 0.22), betatron (P = 0.26), 
and cyclotron (P = 0.19) orbits initiated in the z — 
subspace at the same radial position, r/R — 0.8131, and 
the same energy, 1000 eV, for a simulation time of t = 
10^. Heating of figure-8 orbits occurs at lower Br than 
for cyclotron or betatron orbits. At Br ~ 20 G, FFTs of 
energy for figure-8 ions initiated at higher energy, E > 
O.OSEn ^ 30 keV, in the z = subspace show sharp peaks 
and little further gain in energy. The regular interaction 
with RMFq of these higher energy figure-8 orbits may be 
understood by the greater separation between resonances 
in phase-space [51 12] with increasing energy (or uj). 

We have calculated the PMPo-induced energy gain 
of a figure-8 orbit in a single half period of its motion 
as a first-order correction to the one-dimensional mo- 
tion along r. For P < 0.25, the shape of the effective 
potential, V{r), is a double well[TT]. corresponding to 
cyclotron orbits (in either well) inside the phase-space 
separatrix and to figure-8 orbits (moving across both 
wells) outside the phase-space separatrix. The figure-8 
orbit is approximated by motion in a symmetric double 
well: p = = pq + aicos[u! {t — to)] + a2Cos[3a; it — to)], 
with axjai ~ 10. The amplitudes of oscillation, oi and 
02, are determined by the total energy. The energy 
change from the interaction with the field is given by 
dH/dt = q£ • V = q{£rVr +£^v^), with c£ = —dA/dt. 
After some algebra and integrating over a single half os- 
cillation we get the energy gain from the interaction with 
the £r and components of the RMF. The biggest en- 
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ergy change occurs for a resonance between the RMFq 
and the Fourier components of the orbital motion, i.e., 
lor/uj = s, where s is an integer. The Er- and 5^-induced 
radial and azimuthal portions of the energy change (Ai?) 
in a single half oscillation for an s resonance are: 
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relative fluctuations in energy during an oscillation 
are of order, maxAEodd/ H ^ O (IO^sBr/Bo) , and 
maxAEeven/H ~ O {IQs^Br/ Ba) ■ These predict signif- 
icant energy gain for figure-8 orbits over a single oscilla- 
tion, even for a relatively low amplitude RMF, BR/Ba ^ 
10~^. The energy gain for s — even resonances has an 
dependence while s — odd energy gain has a linear de- 
pendence on s. Resonances with an odd value of s show 
better heating than s — even resonances, especially at 
lower values of s, where the ion energy is higher. Thus, 
the heating observed for figure-8 orbits at higher energies 
results primarily from an overlap of odd-s resonances. 

Using the condition for exponential separation of 
trajectories [12], we now determine the threshold for the 
ergodicity of ion trajectories, essential to convert en- 
ergy gain to stochastic heating. The change in en- 
ergy over an oscillation is used to map the dynamics: 
Ej+i = Ej + AE{tj); tj+i = tj + zn^-^) where 
is the time of the start of successive ion oscillations at 
P = Pmax = Po + ai and AE{tj) is AEritj) + AE^{tj), 
with the substitutions 5* ^ = — ujRtj, and 

AE — > AE{tj). The dynamics will be chaotic if ex- 
ponential separation of trajectories, i.e., K > 1, oc- 
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In dimensionless vari- 



ables, E = {m/b^R^)E and uj — muj/b — 2uj/uici, where 
b = qBa/2c. K for odd and even resonances are: 
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Based on these, increasing the axial wavenumber, fc, of 
the RMFo should lower the chaos threshold for s — odd 



resonances while raising the threshold for s — even reso- 
nances. This is borne out by numerical simulation. 

Fig. 3c) shows dui{E)/dE vs. E for figure-8 orbits hav- 
ing P = 0.15. At energies very close to the separatrix, 
duj{E)/dE grows as {E — Eh)~^^^, where Eh is the en- 
ergy at the phase-space separatrix. The large growth of 
duj{E) / dE near the separatrix corresponds to a large in- 
crease in non-linearity of the figure-8 orbital frequency. 

The greatest rate of stochastic heating is expected to 
occur for lower energy figure-8 orbits where the values of 
s and dCj{E)/dE are higher. As Br is increased, the 
stochastic region above the phase-space separatix will 
broaden. Close to the separatrix, even very low ampli- 
tudes of Br should produce chaotic orbits. Eqns. 3) and 
4), combined with Fig. 3c), can be used to estimate the 
relative amplitude of Br needed to produce stochastic- 
ity and heating. For example, s = 3 resonance occurs at 
E w .0185, corresponding to dCj{E)/dE ss 25, see Fig. 
3c). Using Eqn. 3) and kR ~ 1, chaotic trajectories are 
to be expected for all s — odd resonances with s > 3 and 
BR/Ba > 5 • 10^"*, for ~ 1, the assumption used in the 
derivation. These findings approximately agree with the 
numerical findings. Changing the value of a P changes 
the scale of E, but does not have a substantial effect on 
the value of duj{E)/dE at different resonances. Thus P 
determines the energy range over which figure-8 orbits 
get heated, with greater energy range for lower values 
of P, while not affecting the approximate structure of 
phase space. In Fig. 3c), all s > 4 resonances are located 
to the left of E 0.013, hence occur over the interval 
6E ^ 0.003. This leads to much greater chaos closer to 
the phase-space separatrix where the closely spaced res- 
onances overlap. Thus, lower-energy figure-8 orbits are 
more chaotic and much better heated by the RMF^ than 
the higher energy ones. Fig. 3d) shows this effect for two 
values of P: 0.17 and 0.2. Figure-8 orbits are not further 
heated once their energy reaches (or initially exceeds) the 
curved line appropriate for each P value. The simplifi- 
cations on which Eqns. 3) and 4) are based become less 
accurate at BR/Ba > 0.001. 

Among the clearest differences between these results 
for the FRC and those reported for the tokamak are: 1) 
The non-linearities for the FRC arise from the double po- 
tential well and field gradient and their direct effects on 
the particle orbit. Those in the tokamak arise from trap- 
ping in the wave field - hence require a stronger wave field 
- and resonance between the cyclotron motion and the 
wave field, resulting in a single large resonance and large 
first-order islands. In contrast, close spacing in phase 
space between resonances of a figure-8 orbit leads to an 
overlap between resonances and the observed stochastic 
heating for figure-8 orbits in PRCs. The importance of 
the time-varying field, in the tokamak analysis is to 
create a small nonlinearity in this 2-D system (of the or- 
der of £/B) which leads to resonances between the two 
degrees-of-freedom, not between the £ field and the ion 
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trajectory. 2) In the tokamak, heating occurs at 10 x 
higher values of 17 (over 20 vs 1 in the FRC) and lower 
values of uJa/kvi^thermai 1 vs 10). 3) The threshold 
for heating in the FRC is lower by the factor sdui{E)/dE, 
through which the effect of the PRC's double effective- 
potential well is clear. 

In summary, the energy gain in an orbital period due 
to RMFo was calculated for a figure-8 orbit in an FRC. 
Resonances of w/f with w produce significant energy gain. 
Odd-s resonances more effectively heat for high energy 
(lower s) figure-8 orbits. The energy gain in a oscillation 
was used to map the dynamics and a criterion for the ex- 
ponential separation of trajectories was used to find the 
threshold for chaotic orbits. -ftT, the measure of the rate 
of trajectory separation, increases with Bn- At higher 
energies, the orbits are less chaotic due to both a lower 
value of s and, more importantly, to a decreased nonlin- 
earity reducing dCj{E)/dE. 
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FIG. 1. a) Shape of PRC magnetic field in z-r plane, k = 4.1. 
b) Shapes of typical cyclotron, betatron and figure-8 orbits in 
z = plane. 
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FIG. 2. Results of numerical simulation of ion trajectories 
in an FRC with RMFo- a) Figure-8 ion energy versus time 
for two values of Bh. b) FFTs of ion energy for the two 
cases shown in a) and also for a betatron orbit, c) Separation 
between resonances, A/, versus normalized initial energy for 
P = 0.1, 0.2 and 0.3, at low Br. 




FIG. 3. a) Maximum energy attained versus time for four 
values of Br. b) Maximum energy attained by r = 10* versus 
Br by three types of orbits, c) dw{E)/dE versus E for figure- 
8 orbits. P — 0.15. d) Curves: Threshold E for saturation 
of stochastic heating versus BR/Ba, based on Eqns 3) and 
4): P = 0.17 (upper curve) and P = 0.2 (lower curve); Data 
points, dots and stars, from RMF-code. 



